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Abstract. Large-scale and gap error free calculations of the electronic structure of vacancies 
in 4H-SiC have been carried out using a hybrid density functional (HSE06) and an accurate 
charge correction scheme. Based on the results the carbon vacancy is proposed to be responsible 
for the Z1/2 and EH 6 / 7 DLTS centers. 

Introduction 

The potentials of SiC are compromised by the presence of persistent deep-level defects. Two of 
the most prominent ones are the Z 1 / 2 [1] and the EH 6 / 7 centers [2]. The origin of these centers 
is not well-established but they are both intrinsic-related defects of fundamental nature which 
appear also in as-grown material. Possible candidates are carbon vacancy related defects and 
small carbon clusters at the silicon site. DLTS studies on electron-irradiated 4H-SiC samples 
found that irradiation greatly increases the defect concentrations and the centers persist in 
abundant quantity up to high temperatures. These defects strongly influence charge carrier 
lifetimes which are crucial for bipolar SiC device performance. In particular, due to its negative- 
ly property, the Z 1 / 2 defect is a very efficient electron trap since it can capture two electrons. 

Previous ab initio calculations proposed several microscopic models for the Zi/ 2 center. 
Eberlein et al. suggested a nitrogen-carbon interstitial model. [3] However, later it was shown 
experimentally that this center does not contain nitrogen. [4] Zywietz et al. [5] and Bockstedte et 
al. [6] reported negative-?/ property for the carbon vacancy (Vrj) which might be the Zi/ 2 defect. 
Gali et al. [7] studied carbon-interstitials and proposed a 3-carbon structure showing negative-C/ 
behavior as a candidate for the Zi/ 2 center. Usually, the Z]/ 2 and EH 6 / 7 centers appear together, 
this holds especially for as-grown samples and also after irradiation with low energy (<400 
keV) electrons [8] while after irradiation with MeV electrons or ions the concentrations of the 
two defects deviate from a one-to-one correspondence [9]. In the latter case, the EH 6 / 7 defect 
level is found to be broad with overlapping contributions while in the former case one single 
contribution dominates [8, 10]. Based on these results, Danno et al. [8] suggested that the 
Zi/2 an d EH 6 / 7 levels stem from the same microscopic origin, namely the carbon vacancy. 

Former ab initio simulations on intrinsic defects in SiC were done by standard DFT in 
relatively small supercells with ad-hoc charge corrections, except the study done by Gali et al. 

'] and Bockstedte et al. [6] who utilized a non-selfconsistent method to correct the band gap 
error. Even so, the lack of highly accurate calculations prevented unambiguous identification of 
Zi/2 and EH 6 / 7 defects so far. In this study we have investigated the fundamental monovacancies 
in 4H-SiC and focused on accurate calculation of their occupation levels. 

Method 

It is well-known that standard DFT suffers from band gap error which affects energy position 
of defect levels as well. A highly accurate study needs large supercells and a methodology which 



removes the band gap error as well as corrects charged supercells. We carried out large-scale ab 
initio calculations of silicon and carbon vacancies in 4H-SiC, where the LDA band gap error 
was corrected employing the Heyd-Scuseria-Ernzerhof hybrid functional (HSE06). The HSE06 
functional has been proven to successfully reproduce energy levels for a wide range of defects 
in group-IV semiconductors [11], and the calculated band gap of 4H-SiC of 3.16 eV compares 
well with the experimental value of 3.27 eV. 

First, geometries were optimized by standard DFT in a 576-atom supercell, then the ge- 
ometries were recalculated with the hybrid functional without further relaxation. Calculations 
were performed using the plane-wave projector augmented-wave (PAW) method with Perdew- 
Burke-Ernzerhof (PBE) GGA functional and HSE06 hybrid functional as implemented in the 
VASP code. [12] The plane- wave cutoff for the wavefunction expansion (and for that of the 
charge density) was set to 420 eV (840 eV). For k-point sampling only the T-point was used. 
Besides the band gap error, charged supercells imply another challenge, and we utilized the 
correction scheme suggested by Lany et al. [13] The estimated overall error of the method used 
is around 0.1 eV. 

Results 

The acceptor levels of the carbon vacancy are found to show negative-?/ and close to negative-?/ 
property at the k- and the fa-site, respectively. The difference between the (0/1-) and (0/2-) 
levels remains below 0.1 eV. We did not find any negative-f/ behavior for the donor levels of 
the Vc- The difference between the (2+/1+) and the (1+/0) levels is around 0.1-0.2 eV. The 
positions of the acceptor and donor levels of Vc occure around Ey+(2.7-2.8) eV and Ey+(1.6- 
1.8) eV, respectively, which can explain the levels of the EHg/7 and Z1/2 centers (see Table 1) 
(Ey denotes the valence band maximum). The calculated formation energy of the neutral carbon 
vacancy is around 4.4 and 5.0 eV under Si-rich and C-rich condition, respectively. The cubic 
site is slightly more favorable than the hexagonal one. We also estimated the carbon vacancy 
concentration under typical growth conditions assuming that high temperature CVD growth 
occurs close to quasiequilibrium. Under such conditions, with a growth temperature of around 
1600-1700 °C and the n-doping concentration of around 3-5 xlO 15 cm -3 , the approximated 
carbon vacancy concentration becomes ~2xl0 12 cm -3 which is in the common concentration 
range of Z1/2 and EH6/7 centers in as-grown samples. [14] 

In contrast to standard DFT results, we did not find any donor levels for the silicon vacancy 
(Vsi). Further, previous results have suggested that Vsi can exist even in the 4- charge state [5], 
while our results show that the 3- charge state can only exist in highly n-type material and the 
(3-/4-) level is above the conduction band minimum. The calculated formation energy of the 
neutral silicon vacancy is around 8.0 and 7.5 eV under Si-rich and C-rich condition, respectively. 
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Table 1: Calculated (HSE06) and corrected occupation levels of the vacancies. Values are in eV 
relative to the valence band maximum (Ey). 




Fig. 1: Relative formation energy diagrams of the carbon vacancies (left column) and the silicon 
vacancies (right column). Occupation levels are at the intersection of the lines. The Fermi energy 
(Ep) is relative to E v . 



Hence, for a wide range of doping concentrations the carbon vacancy is much more abundant 
than the silicon vacancy. Only under extreme C-rich conditions and high n-doping levels, the 
formation energy of V|f can approach that of Vc (see Fig. 1). 

Discussion 

The origin of the Z1/2 and EH6/7 centers has been disputed in the literature for a long time. The 
former is a negative-?/ center with a sharp DLTS peak at around Ec-(0.6-0.7) eV, while the 
latter is broader peak with one or more contributions at around Eq- (1.5- 1.6) eV (Eq denotes 
the conduction band minimum). One of the possible microscopic origins for these centers is 
the carbon vacancy, and indeed, our calculated occupation levels of Vc at Ec-(0.5-0.6) eV 
and at Ec-(l-5-1.7) eV are close to the experimentally observed level positions of the Z1/2 and 
EH 6 / 7 defects, respectively. Although the negative-?/ property of the acceptor level remains only 
in the case of the /c-site, the occupation levels of the h-site are very close to each other, within 
0.1 eV, and their negative-?/ property disappears only after applying the charge correction. The 
donor levels of Vc are also close to each other, within 0.2 eV, with no negative-?/ behavior. 
This is consistent with the experimental findings that the EH 6 / 7 peak is broad with one or more 
contributions. The two close donor levels of Vc in the midgap can also explain the appearance 
of the EI5/6 EPR signals, confirming the results by Umeda et al. [15] Nevertheless, it has to be 
noted that for an unambiguous identification of Z 1/2, EH 6 / 7 and EI 5 / 6 centers other microscopic 
candidates have to be investigated as well. 

According to its formation energy, Vsi is a much less abundant defect than Vc, and it 
exhibits only acceptor levels (without negative-?/) in the band gap of 4H-SiC. Its occupation 



levels may explain the acceptor level of the Vsi proposed by Son et al. [16] in semi-insulating 
4H-SiC, with an activation energy around 0.8-0.9 eV, that is close to the calculated (1-/2-) level 
of the Vs; defect at Ec-(0.7-0.8) eV. Furthermore, the calculated properties of Vsi (formation 
energy and level positions) agree closely with those observed experimentally for the so-called 
S-centers [17] and further work is in progress to substantiate this identification experimentally. 

In summary, with a highly accurate ab initio calculations the positions of the occupation 
levels of the carbon and silicon vacancies in 4H-SiC have been estimated. The results suggest 
that the Z1/2 and EH6/7 centers originate from the carbon vacancy. 
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